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FOREWORD 


Several  computer  algorithms  were  recently  developed  at  the  Naval  Surface  Weapons  Center 
(NSWC),  accounting  for  the  gravitational  action  that  the  tidal  dislocations  of  the  earth’s  masses 
exert  upon  artificial  satellites.  The  necessary  physics  background  and  the  computer  routines 
themselves  were  detailed  in  five  technical  reports  between  1976  and  1979.  Included  in  that  body  of 
work  were  a  theory  of  and  a  computer  program  formulation  for  the  perturbing  acceleration  caused 
by  the  ocean  tide.  This  was  based  entirely  on  the  semidiurnal  (M2)  tide  constituent.  Also,  the 
corresponding  solid  earth  tide  term  in  the  equations  of  motion  ignored  ocean  loading. 

Since  then,  an  effort  was  made  to  re-formulate  the  ocean  tide,  introducing  the  Si,  N2,  K2,  Ku 
O ,,  P„  Qi,  Mf,  Mm,  S,„  spectral  constituents.  In  addition,  the  effect  on  the  solid  earth  tide  of 
bottom  loading  by  the  ocean  tide  was  included,  for  convenience,  in  the  perturbing  acceleration  due 
to  the  ocean  tide,  removing  the  need  for  an  extensive  modification  of  the  solid  earth  tide  term. 

The  work  documented  here  was  done  in  the  Space  and  Surface  Systems  Division  in  support  of 
computer  program  development  for  satellite  geodesy.  Joseph  M.  Futcher,  Jr.  of  the  Physical 
Sciences  Software  Branch  is  Programmer.  The  author  is  further  obliged  to  Robert  A.  Manrique  of 
the  SLBM  Research  and  Analysis  Division  for  contributing  computer  results  that  verified  the 
author’s  checkout  calculations. 

This  report  was  reviewed  by  R.  L.  Kulp,  Head,  Space  and  Ocean  Geodesy  Branch;  D.  R.  Brown, 
Jr.,  Head,  Space  and  Surface  Systems  Division;  and  R.  J.  Anderle,  Research  Associate  of  the 
Strategic  Systems  Department. 


R.  T.  RYLAND,  JR.,  Head 
Strategic  Systems  Department 
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INTRODUCTION 


Work  was  recently  undertaken  to  establish  theories  and  computer  algorithms  for  the 
gravitational  action  by  which  the  tidal  redistributions  of  the  earth’s  masses  perturb  satellite  orbits. 
From  this  effort  resulted  five  technical  reports1-5  of  which  the  first  four  document  the  physics 
background  and  details  of  the  tide  models.  The  fifth  specifies  the  computer  routines  for  the  various 
perturbing  accelerations  resulting  from  these  models,  thus  permitting  the  tide  effects  to  be  included 
into  the  equations  of  motion  of  the  CELEST  and  TERRA  computer  programs. 

Contained  among  the  tidal  perturbing  terms  was  an  algorithm  (Chapter  3,  Pages  17  through 
31,  of  Reference  5)  for  the  ocean  tide.  This  started  from  an  existing  table  for  the  global  tide 
amplitudes  and  phase  angles6-8,  postulating  that  each  tidal  height  value  resembles  a  gravitating 
point  mass  that  by  its  presence  perturbs  the  satellite  motion.  Only  the  semidiurnal  (A/2)  lunar  tide 
was  considered.  But  this  algorithm  was  regarded  as  setting  the  pattern  for  introducing  other 
components  of  the  fluid  tide  later.  The  restriction  has  since  been  removed  because  amplitude  and 
phase  angle  tables  have  become  available  for  the  Slt  N2,  Klt  Kt,  O,,  Pu  Qu  Mf,  Mm  and  S,„ 
constituents  of  the  ocean  tide.  It  then  became  possible  to  revise  the  existing  computer  program  for 
the  perturbing  acceleration  due  to  the  A/2  tide  by  formulating  the  perturbing  terms  for  the  just 
mentioned  additional  members  of  the  sea  tide  spectrum. 

Part  of  the  1979  collection  of  tide  algorithms5  was  a  computer  program  formulating  the 
Newtonian  attraction  caused  by  the  tide  bulge  of  the  solid  earth.  This  included  both  the  lunar  and 
the  solar  tide  components.  The  tidal  properties  were  assumed  to  be  functions  of  latitude. 
Accordingly,  Love  coefficients  were  featured  that  are  zonal  expansions  of  latitude.  An  allowance 
was  made  for  tide  lag  manifesting  itself  as  a  time  delay  in  the  response  of  the  tidal  mass 
redistribution  to  the  tidal  stress  field.  Absent  was  however  the  effect  of  ocean  tide  loading  on  the 
tide  bulge.  It  appeared  soon  desirable  to  account  for  ocean  tide  loading.  To  avoid  having  to  burden 
the  algorithm  for  the  land  tide  with  the  very  numerous  parameters  specifying  the  ocean  tide,  it  was 
decided  (and  found  possible)  to  leave  the  former  unchanged  and  to  reflect  the  necessary  alterations 
in  the  computer  routine  for  the  latter.  This  happened  to  be  not  only  the  most  economic  but  also, 
from  the  viewpoint  of  the  tide  physics  involved,  the  most  natural  way  of  doing  things.  The  task 
turned  out  to  consist  merely  of  a  trivial  modification  of  the  tidal  point  mass  in  the  ocean  tide  model. 


REVISED  COMPUTER  ALGORITHM  FOR  THE  PERTURBING  ACCELERATION 
CAUSED  BY  THE  OCEAN  TIDE  INCLUDING  THE  OCEAN  LOADING  EFFECT 
ON  THE  SOLID  EARTH  TIDE  ROUTINE 


INPUT  DATA 

Concerned  are  the  ocean  tide  modes  Mu  S2,  Nt,  Ku  Kit  O,,  P„  Qlt  Mf,  Mm  and  S,„.  For  each 
of  these  tides,  there  exists  a  computer  tape  that  lists  the  NP  values  of  the  tidal  amplitude  and  the 
tidal  phase  6if. 


T 

A IP  Number  of  grid  points.  This  may  be  expected  to  be  a  five-digit  integer. 

Z.j  Tidal  amplitude  on  area  element  ij.  If  available  in  meters,  use  directly.  | 

Otherwise,  convert  to  meters.  I 

6,j  Tidal  phase  angle  on  area  element  ij,  in  degrees. 

Custodian  of  the  tapes  just  mentioned  is  Ling  Szeto  of  the  Physical  Sciences  Software  Branch.  Also 
required  are 


R  “Radius  of  the  Earth,’’  i.e.,  the  semimajor  axis  of  a  suitable  reference 
ellipsoid,  in  kilometers.  Trial  value:  R  =  6378.145  km. 

e2  Square  of  eccentricity  of  reference  ellipsoid.  Trial  value:  c2  =  0.00669342. 
In  case  it  is  desired  to  start  from  the  ellipsoid  flattening,/,  find  e2  from  t2  = 
(2  -  /)/ 

Me  Gravitational  constant  of  the  earth,  in  kilometers  and  seconds.  Trial  value: 
Me  =  398601  km3  sec'1. 

G  Newton’s  constant,  in  km,  kg  and  sec. 

Trial  value:  G  =  6.6732E-20  km3  kg' '  sec'2. 

6  Density  of  water,  in  kg  and  km. 

Trial  value:  q=  1,  F.+  \2kgkm'3. 

q*  Density  of  the  ocean  bottom,  in  kg  and  km. 

Trial  value:  q*  =  3. E  +  12  kg  km'3. 

o  Rate  of  mean  mean  longitude  of  moon. 

a  =  J80  1 .40519E-04  deg  sec'. 
n 

The  following  parameters  J,  n  and  t*  speci  .  the  time  instant  at  which  the  perturbing  acceleration 
due  to  the  tide  is  to  be  found.  Normally,  this  is  the  time  instant  associated  with  the  current  time  line 
of  the  orbit  integration. 

J  Number  of  the  calendar  year. 

n  Number  of  the  day  within  the  year. 

t *  Mean  solar  time  at  Greenwich  (GMT,  UTC),  in  seconds. 
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xt  )  Inertial,  Cartesian  components  of  satellite  position  vector  in  km,  associated 
Xi  >■  with  the  just  specified  time  instant.  For  the  precise  definition  of  the  reference 
x3  J  frame,  see  TERRA  and  CELEST  documentation9- 

( ABCD)j  „ Transformation  matrix  that  manages  the  transition  from  the  inertial  frame 
of  the  satellite  equations  of  motion  (“Basic  Inertial  System,”  which  is 
associated  with  either  1950.0  or  with  the  beginning  of  the  day  UT  of  the 
trajectory  epoch)  to  the  earth-fixed  reference  frame  corresponding  to  the 
time  instant  J,  n,  t*.  This  transformation  is  a  computer  routine  that  is 
already  part  of  TERRA  and  CELEST.  For  details,  see  the  TERRA  and 
CELEST  documentation9. 

NMAX  Limit  on  the  number  of  terms  in  the  expansion  of  the  tide  potential. 


COMPUTER  ROUTINE  FOR  THE  EXPANSION  COEFFICIENTS 
OF  THE  TIDE  POTENTIAL 

The  time-independent  constituents  of  the  expansion  coefficients  (aF„m,  PF„m,  aHnm,  /*//„„)  of 
the  tide  potential  will  now  be  calculated  from  the  amplitude  and  phase  angle  computer  tapes, 
separately  for  each  of  the  eleven  ocean  tide  modes. 

Once  established,  the  coefficients  aFnm,  PF„„,  aHnm,  PHnm  will  remain  unchanged  as  long  as  the 
amplitude  and  phase  tape  for  the  particular  tide  mode  remains  valid.  Thus,  the  computer  routine 
under  discussion  will  be  exercised  quite  infrequently,  namely  just  once  each  time  a  new  edition  of 
the  just  mentioned  amplitude  and  phase  tape  becomes  available.  That  is  expected  to  occur  at 
intervals  of  several  years,  during  which  the  present  routine  will  remain  dormant. 

For  each  ocean  tide  mode  execute  the  following  algorithm.  For  the  geometry  associated  with 
the  index  ij,  see  Figures  1,  2  and  3. 

a  a  =  IQ'3  (e  -  0.0667c  *)  G  AS*  £*  cos  d*  (01 ) 


pij  =  lO'3  (g  -  0.0667C  * )  G  AS*  £*•  sin  d*  (02) 


=  area  of  the  non-polar  surface  area  element  for 
j  —  2,  3,  4,  •  •  ,  jMAx  <  180 

(as*)'™'4*  =  y  {~m)  R1  (04) 

=  area  of  the  polar  surface  area  element  (neighboring  the 
North  Pole) 


Figure  2.  Position  of  Point  Mass  at  the  Geometrical  Center  of  the 
Surface  Area  Element  Associated  With  ij 


y1  *  INDICATES  THE  POINT  WHERE  GREENWICH  MERIDIAN  AND  EQUATOR  INTERSECT 
y2  *  LOCATED  IN  THE  EQUATORIAL  PLANE 

Figure  3.  Position  of  the  Point  Mass  M,}  in  the 
Earth-Fixed  Cartesian  Coordinate  Frame 
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02) 


d*  = 


I  =  k 
t  *  k 


Calculate  the^m  and  h"m  as  follows,  noting  that 


(13) 


and 


and  that 


hvno  =  0  for  all  values  of  n 


(14) 


Km  •Am  =  0 for any n<m-  (15) 

Obtain,  separately  for  each  v,  the  required pnm  and  hvnm  from  the  following  recurrence  relations: 

To  advance  in  n,  evaluate 

fi  +  u*  =  n-  m  +  i  K2w  +  l)sin Qvfv„,m  -  («  +  m)pvA„J  (16) 

=  ~ — ^rr t[(2«+  l)s>n9v^.m  -  (n  +  m)pvh\_,  J  (17) 

n  —  m  +  1 


To  advance  in  m,  use 

=  (2*  +  1 )  A> (cos  9V  cos  KA„.„  -  cos  9V  sin  Av  hv„  n] 
K^.n.x  =  (2w  +  1  )Pv  [cos  9V  cos  A„ hvn„  +  cos  0„  sin  Xvp„ J 
Start  the  recurrences  from 
fv  =  — L 

J  0.0  ^ 


A,  =  * 


sin  8V 


(18) 

(19) 


(20) 

(21) 


<*  =  <0  =  0  (22) 
and  terminate  the  procedure  at  n  =  NMAX. 

Store  the  resulting  0Fnm,  ®Fnm,  aH„m  and  PHnm.  They  will  be  programmed  as  constant 
parameters  into  the  computer  routine  for  the  perturbing  acceleration  and  are  expected  to  serve  for 
all  subsequent  computer  runs  of  that  routine,  until  updated.  Note  that  the  °F„„,  PFnm,  °H„m  and 
pH„m  are  dimensionless  quantities. 
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CALCULATION  OF  THE  ASTRONOMICAL  ARGUMENTS 
FOR  THE  VARIOUS  TIDE  COMPONENTS 


Evaluate  now,  using  double  precision  throughout. 


6(n,m)  = 


1  n  =  m 
0  ^  n  +  m 


Nt  =  n  <5(7,  1975)  +  (365  +  n)  <5(7,  1976) 

+  (  731  +  n)  <5(7,  1977)  +  (1096  +  n)  6(J,  1978) 

+  (1461  +  n)  <5(7,  1979)  +  (1826  +  n)  <5(7,  1980) 

+  (2192  +  n)  <5(7,  1981)  +  (2557  +  n)  <5(7,  1982) 

+  (2922  +  n)  <5(7,  1983)  +  (3287  +  n)  <5(7,  1984) 

+  (3653  +  n)  <5(7,  1985)  +  (4018  +  n)  <5(7,  1986) 


A T  =  (5.28E-04)  +  (3.56£-08)  Nx 


d0  =  27392.5  +  N*  +  AT 


0  36525 

h0  =  279.69668 

+  36000.7689304850  r0 
+  0.000303  Tl 


So  =  270.434358 

+  481267.88314137  T0 
-  0.001133  Tl 
+  0.0000019  Tl 


(101) 


(102) 

(103) 

(104) 

(105) 

(106) 

(107) 


p0  =  334.329653 

+  4069.0340329575  T0 

-  0.010325  Tl 

-  0.000012  Tl  (108) 

For  each  tide  mode,  find  now  the  Astronomical  Argument  x  fr°m 

Tide  Mode  x  *n  deg  = 

Mi  2  (ho  —  So) 

Si  0 


Ni 


2  h0  -  3  So  +  Po 
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(109) 

(110) 
(111) 


Tide  Mode 


X  in  deg  = 


Ki 

2  h0 

h0  +  90 

o, 

ho  —  2  5o  90 

P, 

h0  -  90 

e. 

h0  -  3  s0  +  Po  -  90 

M, 

2  So 

So  —  Po 

s„ 

2  h0 

(112) 

(113) 

(114) 

(115) 

(116) 

(117) 

(118) 

(119) 


Note  that  the  day  number  “n,  ”  and  thus  x,  niust  be  updated  whenever  the  time  argument  t*  runs 
into  the  next  day. 


COMPUTER  ROUTINE  FOR  THE  PERTURBING  ACCELERATION 

For  each  tide  mode,  separately,  use  the  set  of  time-independent  tide  potential  coefficients  and 
the  astronomical  tide  argument  to  evaluate  the  following  algorithm. 

First,  calculate  the  earth-fixed  Cartesian  satellite  coordinates  (y„  y2,  y3)  from  the  corresponding 
inertial  coordinates  (*„  x2,  jc3),  the  latter  being  associated  with  the  time  line  at  t*  seconds  from 
epoch. 


(201) 


(202a) 


(202b) 


as  convenient. 

Calculate  now  the  time-dependent  expansion  coefficients  for  the  tide  potential  (note:  n,  m  < 
NMAX)  from 


Fnm  =  aFnmcos(ot*  +  x)  +  pF„m  sin  (of  +  x) 


(203) 


Hnm  =  °H»m  cos  (of  +  x)  +  pH„„  sin  (of  +  x)  (204) 

The  argument  (of  +  x)  should  be  evaluated  in  double  precision.  The  double  precision  sin  and  cos 
functions  should  be  applied.  Then  calculate  Fnm  and  Hnm  and  subsequently  revert  to  single  precision. 

Note  also  that  F„m  and  H„m  are  linear  functions  of  the  two  trigonometric  terms.  To  evaluate  these 
trigonometric  terms,  let  /*and  /,* ,  be  the  values  of  Universal  time  for  which  subsequent  integration 
steps  are  to  be  performed.  To  save  computer  time,  update  the  trigonometric  time  factors  as  follows: 


cosfo/,*,  +  x)  =  cos[(o/*+  x)  +  oA/]  =  cos  (ot*+  x)cosoAf  -  sin(of,*+  x)sinoAf  (205) 

sin  (of*  i  +  x)  =  sin  [(of  *+  x)  +  oA/]  =  sin(o/*+  x)cosoAf  +  cos  (of,*+  x)sinoAf  (206) 

A/  =  t*t  ,  -  t*  (207) 

Now  calculate  the  Cartesian  components  of  the  perturbing  acceleration  in  the  earth-fixed  frame: 

3«t> 


where 


dy. 


NMAX  n 

I  s 


n=0  m= 0 


NMAX 


II 

z. 

n  ~0 

z. 

m  =  0 

[F„m 

v  w  nm 

dy. 

+  Hnm 

”  r  nm 

dy. 

3* 

NMAX 

y 

n 

y 

If 

dUnm 

+  Hnm 

dVnm 

dy. 

n  =  0 

m-0 

l  *  nm 

dy. 

dy , 

dU*m 

dy, 

du„m 

dy, 

dUnm 

dy, 

dv„m 

dy, 


1  / _ L  a  y  _ L  V  \ 

JR  t  2  2  «♦!  ,m*l  J 


-1 


(n  -  m  +  \)Un. 


_L  ( —  A  V  _ -V  \ 

^  V  2  mn  r  n  *  *  * m  ■  *  2  n  *  *  ’ m  * 1 J 


(208) 


(209) 


(210) 


(211) 

(212) 

(213) 

(214) 
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ay 1  R  \2 


«*  mn 


1.  m  t 


2 


fl  +  I,  m  +  / 


A  1/  _  1 

— ^  =-5L(«  -  m  +  1)  Kntl,m 

9^3  /? 

(216) 

and 

Am„  -  (n  -  m  +  l)(n  -  m  +  2) 

(217) 

n  _  («  -  D!  _  Uml 

'  (n  +  1)!  n{n  +  1) 

(218) 

,/  _  .  -  D!  p  _  _  K„ 

(n  +  l)f  n(n  +  1) 

(219) 

Obtain  the  values  of  the  individual  eigenfunctions  from  the  following  recursive  relationships: 

R 

P-~ 

(220) 

V3 

sin  yj  =  -Aj- 

(221) 

V„ o  =  0  for  all  values  of  n 

(222) 

Unm  --  V„m  =  0  for  all «  <  m 

(223) 

To  advance  in  n,  evaluate 

c/- .  1. »  =  - - SrXT  +  1 )  sin  t+»  U„m  -  {n  +  m)pU„.tm\ 

n  —  m  +  i 

(224) 

=  - ~r~ r  K2"  +  1 ) sin  1*3  Vnm  -  (n  +  m)/? 

rt  —  W  +  1 

(225) 

To  advance  in  m,  use 

V„  (2/t  +  l)p  l/„„  - 

(226) 

K. .  .  ,  =  (2n  +  \)P  (-^f-  K,„  +  «/„) 

(227) 
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Start  from 


t/o.0  = 

r 

(228) 

(/,„  = 

Ry> 

- 

(229) 

rJ 

K,o  = 

V,.o  =  o 

(230) 

The  Cartesian  components  of  the  perturbing  acceleration,  in  the  earth-fixed  reference  frame,  are 
now 


T  -JL 
"  dyt 


(231) 


T  = 


1  ?2 


dy2 


T  = 


a<j> 


dy, 


(232) 

(233) 


Finally,  rotate  the  perturbing  acceleration  vector  back  into  the  inertial  coordinate  system: 


T 

rx2  )  =  (ABCD)l,.  (  Tn 


(234) 


where  ( ABCD)T  is  the  transpose  of  (A BCD). 


COMPUTER  PROGRAM  OUTPUT 

Add  up  thex,  components,  Tx ,,  of  the  perturbing  accelerations  resulting  from  the  eleven  tide 
modes 


Tx  i  —  (Txi)Mi  +  (Tx,)S2  +  •••••+  (7',i)j,0 

and  also  compute  for  the  x2  and  x3  components 

T,l  =  ( Txl)m  +  •••••••••  +  ( Txl)s,a 

Txj  =  (T,  j)M2  + . .  (Tx])s.a 


(301) 


(302) 

(303) 
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These  are  the  Cartesian  components,  in  inertial  space,  of  the  perturbing  acceleration  caused  by  the 
combined  action  of  all  ocean  tide  modes  considered  plus  the  effect  on  the  solid  earth  tide  of  tidal 
ocea  .  loading  (which,  as  already  mentioned,  had  not  been  included  in  the  formulation  of  the 
perturbing  acceleration  produced  by  the  solid  earth  tide  in  the  previous4  5  tide  routines). 


REFERENCES 

1.  W.  Groeger,  “A  Gravitational  Potential  for  Atmospheric  Earth  Tides  Caused  by  the  Sun”, 
Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL  TR-3485, 
Dahlgren,  VA,  October  1976. 

2.  R.  Manrique  and  W.  Groeger,  “A  Gravitational  Potential  for  Atmospheric  Earth  Tides 
Caused  by  the  Moon”,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical 
Report  NSWC/DL  TR-3638,  Dahlgren,  VA,  November  1976. 

3.  W.  Groeger  and  R.  Manrique,  “ Ocean  Tides  in  the  Equations  of  Motion  of  the  Terra  System 
of  Computer  Programs  for  Satellite  Geodesy”,  Naval  Surface  Weapons  Center,  Dahlgren 
Laboratory  Technical  Report  NSWC/DL  TR-3535,  Dahlgren,  VA,  December  1976. 

4.  W.  Groeger  and  R.  Manrique,  ‘‘Solid  Earth  Tides  in  the  Equations  of  Motion  of  the  Terra 
System  of  Computer  Programs  for  Satellite  Geodesy”,  Naval  Surface  Weapons  Center, 
Dahlgren  Laboratory  Technical  Report  NSWC/DL  TR-3560,  Dahlgren,  VA,  February  1977. 

5.  W.  Groeger,  “Five  Algorithms  for  the  Earth  Tide  Terms  in  the  Force  Model  of  NS WC 
Computer  Programs  for  Satellite  Geodesy”,  Naval  Surface  Weapons  Center,  Dahlgren 
Laboratory  Technical  Report  NSWC/DL  TR-3912,  Dahlgren,  VA,  February  1979. 

6.  E.  W.  Schwiderski  and  L.  T.  Szeto,  “NSWC  Ocean  and  Geocentric  Tide  Data  Tapes  and  Tide 
Computation  Program ”,  Preliminary  Report  (Naval  Surface  Weapons  Center,  Dahlgren 
Laboratory  Technical  Report  manuscript),  Dahlgren,  Virginia,  15  September  1979. 

7.  E.  W.  Schwiderski,  “Global  Ocean  Tides,  Part  I:  A  Detailed  Hydrodynamical  Interpolation 
Model”,  Naval  Surface  Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL 
TR-3866,  Dahlgren,  VA,  September  1978. 

8.  E.  W.  Schwiderski,  “Global  Ocean  Tides,  Part  II:  The  Semidiurnal  Principal  Lunar  Tide 
(M2),  Atlas  of  Tidal  Charts  and  Maps”,  Naval  Surface  Weapons  Center.  Dahlgren  Laboratory 
Technical  Report  NSWC/DL  TR-79-414,  Dahlgren,  VA,  December  1979. 

9.  J.  W.  O’Toole,  “Celest  Computer  Program  for  Computing  Satellite  Orbits”,  Naval  Surface 
Weapons  Center,  Dahlgren  Laboratory  Technical  Report  NSWC/DL  TR-3565,  Dahlgren, 
VA,  October  1976. 


12 


10.  E.  W.  Schwiderski,  "Ocean  Tides,  Par,  !:  Global  Ocean  Tide  Equations",  1.  Marine  Geodesy, 
Volume  3,  1980. 


APPENDIX  A 

NOTES  ON  AUGMENTATION  OF  SOLID  EARTH  TIDE  ALGORITHM 

BY  OCEAN  LOADING 


15 


NOTES  ON  AUGMENTATION  OF  SOLID  EARTH  TIDE  ALGORITHM 

BY  OCEAN  LOADING 


The  outstanding  structural  difference  between  the  above  documented  algorithms  for  the 
expansion  coefficients  of  the  tide  potential  and  perturbing  acceleration  and  the  corresponding 
formulae  and  computer  routines  in  References  3  and  5  is  that  the  former  account  for  eleven  tide 
modes  (including  the  M2  tide)  while  the  latter  are  restricted  to  the  semidiurnal  tide.  Compared  with 
it,  the  modification  that  now  introduces  the  influence  of  hydrostatic  ocean  tide  loading  on  the  solid 
earth  tide  appears  rather  minor.  In  fact,  it  consists  of  just  a  slight  alteration  of  the  equations  for  the 
quantities  a,j  and  fi0  (Equations  01  and  02  in  the  main  part  of  this  report).  The  following  is  the 
theoretical  background  of  this  alteration*). 

For  each  tide  mode,  the  tidal  elevation  of  the  sea  surface  is 
4»  =  L  cos  [(at*  +  x)  - 
=  £„  cos  d„  cos  (at*  +  x) 

+  4„sind„  sin  (at*  +  x)  (401) 

where  the  astronomical  tide  argument,  x.  is  a  linear  function  of  h0,  s0  and  p0 

X  =  fct[h0(To),s0(T0),Po(T0)]  (402) 


Figure  4.  Geometrical  Relationship  Between  Ocean  Tide  and  Sea  Bottom  Tide. 


*  Note  that  except  for  definitions  of  the  quantities  occurring  are  unnecessary  because  they  may  be  readily  found  in  the 
main  part  of  this  report  as  well  as  in  References  3  and  5.  As  before,  units  are  km  for  length,  kg  for  mass  and  sec  for  time 
except  for  the  tidal  amplitude  on  area  element,  i0,  specified  in  cm. 


T0  is  the  fractional  time  in  centuries  at  the  beginning  of  the  (mean  solar)  day  UT.  The  point  mass  m. 
representing  that  portion  of  the  tide  bulge  which  is  located  within  the  particular  surface  area  element 
AS.  (Reference  3,  Page  3)  is 


my  =  (o.  cos  (of*  +  x)  +  0,  sin  (of*  +  x)J 

o 

(403) 

a,  =  qG  AS.  10~35.cosd. 

(404) 

0.  =  pGAS.  10'3£.sind. 

(405) 

Defining  a  “gravitational  constant”  ji,  for  each  point  mass  one  can  also  say 
fj.  =  Gmy  =  a.  cos  (of*  +  x)  +  /?.  sin  (of*  +  x) 


(406) 


If  one  now  recalls3-5  that  the  perturbing  potential  aloft  caused  by  the  ocean  tide  may  be  expressed  as 


* 

NMA  X  n 

=  Z  Z  [FHmUmm  +  HnmVnm] 

n  -  0  m  =  0 

(407) 

Fm 

=  i  p 

nm 

v  =  1 

(408) 

Hnm 

N 

=  Z  H'nm 

V  =s  ] 

(409) 

F\m 

_  (2  Aoj  M.  («  -  "»)!  — fir — PT(sin  0.)  cos  m\y 

H  (n  +  m)l  R’  ' 

(410) 

HL 

=  -Hr-  — |^-P:(sin0.)cosmA. 

M  (n  +  m)\  R" 

(411) 

uHm 

=  — — P^(sin  9)  cos  mX 

(412) 

vnm 

=  — — p^(sin  6)  cos  mX 

(413) 

it  becomes  evident  that  it  is  the  quantities  fi,  which  introduce  time  into  the  expansion  coefficients  Fnm 
and  Hnm  of  the  potential. 

To  realize  how  the  effect  of  ocean  tide  loading  may  be  formulated,  consider  Figure  4.  This 
illustrates  the  interaction,  by  hydrostatic  loading,  of  the  ocean  tide  with  the  solid  earth  tide.  The 
present  approach  is  based  on  the  fortunate  circumstance  that,  as  evidenced  by  Reference  6  (see 
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especially  Equation  8),  Reference  10  and  Reference  8  (see  Equation  3),  the  dip  zEO  of  the  solid  earth 
(ocean  floor)  in  response  to  the  oceanic  tidal  load  is  proportional  to  the  ocean  tide  elevation 

ZST  -  Za  =  Z£°  ~  +  0.0667  |  (414) 

One  must  further  keep  in  mind  that  the  perturbing  acceleration  based  on  the  solid  earth  tide,  as 
previously  formulated4  5,  contains  the  effect  of  that  portion  of  the  oceanic  crust  that  would  in  the 
absence  of  hydrostatic  ocean  tide  loading  be  located  between  the  surfaces  marked  “Solid  Earth 
Tide”  and  “Bottom  Tide”  in  Figure  4.  Because  the  sea  floor  is  now  depressed  by  ZEO ,  we  must 
nullify  the  just  mentioned  part  of  the  solid  tide.  This  may  conveniently  be  done  by  asserting  that  the 
missing  mass  is  a  negative  mass  layer  whose  gravity  acting  on  the  satellite  will  offset  the  unwanted 
masses  present  in  the  existing  version  of  the  solid  earth  tide  algorithm. 

Making  use  of  the  fact  that  the  ocean  depth  is  a  small  fraction  of  the  earth’s  radius,  the 
perturbing  acceleration  due  to  this  negative  mass  layer  (density:  q *,  located  on  the  ocean  surface) 
may  now  be  derived  in  exactly  the  same  manner  in  which  the  perturbation  due  to  the  sea  tide  was 
developed  in  Reference  3.  There  will  result  a  perturbing  potential  identical  to  Equations  403  through 
413  above,  except  that  the  water  density  q  is  replaced  by  -  0.0667  g*  (see  “Input  Data”  in  the  main 
part  of  this  report).  The  combined  gravitational  action  of  the  ocean  tide  and  the  just  defined 
“missing  mass  layer”  is  then  reflected  in  the  main  part  of  this  report,  the  only  difference  between  it 
and  the  case  of  a  pure  ocean  tide  consisting  of  the  term  -0.0667  g*  in  Equations  01, 02, 415  and  416, 

a.  =  (g  -  0.0667  g*)  G  AS,  103  £,  cos  d,  (415) 

0,  =  (g  -  0.0667  g*)  G  AS,  10‘3  £,  sin  d,  (416) 

Otherwise,  the  ocean  tide  algorithm  remains  undisturbed. 
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TRIAL  DATA  FOR  COMPUTER  PROGRAM  CHECKOUT 


The  following  trial  calculations  are  intended  to  facilitate  computer  program  checkout.  They 
were  done  on  an  electronic  calculator.  To  be  feasible,  they  were  restricted  to  the  lunar  semidiurnal 
(M2)  tide  mode.  Also,  for  simplicity,  the  density  of  the  ocean  floor  rock,  e*,  was  assumed  zero. 


INPUT  DATA 

NP  =  9 


Quantities  Associated  with  the  Surface  Area  Elements 


ij 

V 

i, 

m 

<5 

deg 

e 

rad 

A 

rad 

Q 

km 

AS„ 
km 2 

0tV  fiy 

km*  sec~2  km3  sec1 

1,1 

1 

10 

25 

1.562069681 

8.7266463  E-03 

6356.800824 

108.1411251 

6.5403462E-08  3.04981 35E-08 

2,1 

2 

10 

25 

1.562069681 

2.61799388E-02 

6356.800824 

108,1411251 

6.5403462E-08  3.0498 135E-08 

3,1 

3 

10 

25 

1.562069681 

4.36332313E-02 

6356.800824 

108.1411251 

6.5403462E-08  3.0498135E-08 

1,2 

4 

10 

25 

1.544616388 

8.7266463  E-03 

6356.813825 

432.4766612 

2.6156072E-07  1.2196777E-07 

2,2 

5 

20 

30 

1.544616388 

2.61799388E-02 

6356.813825 

432.4766612 

4.9987043E-07  2.8860033E-07 

3,2 

6 

20 

30 

1.544616388 

4.363323 12E-02 

6356.813825 

432.4766612 

4.9987043E-07  2.8860033E-07 

1,3 

7 

10 

25 

1.527163095 

8.7266463  E-03 

6356.839812 

648.550316 

3.9224149E-07  1.8290521E-07 

2,3 

8 

20 

30 

1.527163095 

2.61799388E-02 

6356.839812 

648.550316 

7.496153  E-07  4.327906  E-07 

3,3 

9 

20 

30 

1.527163095 

4.363323 13E-02 

6356.839812 

648.550316 

7.496153  E-07  4.327906  E-07 

R  =  6378.145  km 

£2  =  0.00669342 

Me  =  398601  km3  sec'2 

G  =  6.6732E-20  km3  kg"  sec'2 

e  =  l.£+  12  kg  km'3 

Q*  =  0 

180 

o  - -  1.40519E-04  deg  sec'1  (evaluate  in  double  precision) 

rr 

—  =  .318309886183790671538 

n 

23 


.  ->■ 


J  =  1977 


n  =  202 

t*  =  50000  sec 

xt  =  +  3151.52923  km 

X!  =  +  5458.60875  km 

X3  —  +  3639.07250  km 

/All  A12  A13\ 

( ABCD)jnJ .  =  A21  A22  A23 

\A31  A32  A33/ 

All  =  -  .8405285753 

A12  =  +  .5417623775 

A13  =  +  .2289080162  E-02 

A21  =  -  .5417605355 

A22  =  -  .8405316908 

A23  =  +  .1413662999  E-02 

A31  =  +  .2689913850  E-02 

A32  =  -  .5190827376  E-04 

A33  =  +  .9999963803 

NMAX  =  4 
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COMPUTATION  OF  TIME-INDEPENDENT  POTENTIAL  COEFFICIENTS 


Values  for  the  F  and  h" 

rt  m  n  m 


V 

<*0,0 
km  1 

h0.0 
km  1 

F 

1,0 

km' 1 

h«.o 
km' 1 

I 

1.5731 184E  -04 

0 

I.5783403E-04 

0 

2 

1.5731 184E-04 

0 

1.5783403E-04 

0 

3 

1.5731 184E- 04 

0 

1.5783403E-04 

0 

4 

1.5731  15  IE  -  04 

0 

1 .5778531E  -  04 

0 

5 

1.5731  15  IE -04 

0 

1 .577853 1 E  -  04 

0 

6 

1 .573 1  1 5  IE  -  04 

0 

1 .577853 1 E  -  04 

0 

7 

1.5731087E  -04 

0 

1.5768788E  -04 

0 

8 

1.5731087E  -  04 

0 

1 .5768788E  -  04 

0 

9 

1.5731087E  -  04 

0 

1 .5768788E  -  04 

0 

V 

km-1 

K  i 

I  i  * 

km'1 

F 

2,0 

km'  1 

h2.0 
km'  1 

1 

1 .3773442E  -  06 

1.20 1990 IE  -08 

1.5835193E  -  04 

0 

2 

1.3769246E  -  06 

3.6056041  E -08 

1.5835193E  -04 

0 

3 

1.3760857E  -  06 

6.0081 198E- 08 

1.5835193E  -  04 

0 

4 

4.13 15963E  -  06 

3.6055895E  -  08 

1.5820627E  -  04 

0 

5 

4.1303378E  -  06 

1.0815670E  -07 

1 .5820627E  -  04 

0 

6 

4.127821  IE  -06 

1.8022456E  -  07 

1 .5820627E  -  04 

0 

7 

6.8845393E  -  06 

6.0080464E  -  08 

1 .579 1 5 13E  -  04 

0 

8 

6.8824422E  -  06 

1.8022309E  -  07 

1 .57915 13E  -  04 

0 

9 

6.8782468E  -  06 

3.003 1082E -07 

1. 57915 13E  -04 

0 

V 

F 

‘2,1 

km' 1 

h2.1 

km'1 

F 

2,2 

km'1 

h2.2 

km'1 

I 

4. 1457488E  -  06 

3.61 79402E  -  08 

3.6175267E  -  08 

6.3I44I63E-  10 

2 

4. 1444859E  -  06 

1. 08527 18E  -  07 

3.613 1 193E  -  08 

1.8935556E-09 

3 

4. 1419607E  -  06 

1.8084 19 !  E  -  07 

3.6043099E  -  08 

3.1533625E-09 

4 

1.24321  20E- 05 

I.0849347E  -  07 

3.2550933E  -  07 

5.6817864E  -09 

5 

1 .2428333E  -  05 

3.2544735E  -  07 

3.251 1274E  -  07 

1 .7038437E  -  08 

6 

1 .2420760E  -  05 

5.42302 10E  -07 

3.2432006E  -  07 

2.8374329E  -  08 

7 

2.0703 116E-  05 

1.8067335E  -  07 

9.038143 1 E  -  07 

1.5776137E-08 

8 

2.0696809E  -  05 

S.4196503E  -  07 

9.0271 3 1 5E  -  07 

4.7309 1 91 E -08 

9 

2.0684 198E- 05 

9.0309 16 IE -07 

9.005 1217E- 07 

7.8784606E  -08 

25 


>0  Oo  S'  w  ^  U  IJ 


Values  for  the  F  and  h*  (Continued) 

n  in  n  m 


V 

fV 

*3,0 

km-1 

h3.0 

km" 1 

F 

‘3.1 

km"1 

h3 ,1 
km  1 

1 

1.5886548E  -  04 

0 

8.3188628E  -  06 

7.259761 5E- 08 

2 

1.5886548E-04 

0 

8.3 163287E  -  06 

2. 1777073E  -  07 

3 

1.5886548E  -  04 

0 

8.31  12614E-06 

3.628775  IE -07 

4 

1 .585739 1 E  -  04 

0 

2.493485 IE  -  05 

2. 1 7603 1 5E  -  07 

5 

1.585739 IE  -  04 

0 

2.4927255E  -05 

6.52743 16E-  07 

6 

1. 5857391 E -04 

0 

2.491 2067E  -  05 

1.0876843E  -  06 

7 

1.5799154E  -  04 

0 

4. 1485684E  -  05 

3.6204008E  -  07 

8 

1.5799154E  -  04 

0 

4. 1473047E  -  05 

1 .0860 1 00E  -  06 

9 

1.5799154E  -  04 

0 

4.I447777E  -05 

1  8096490E  -  06 

V 

F 

3,2 

km" 1 

h3,2 
km  1 

F 

*3.3 

km  1 

h3.3 

km"1 

1 

1.8147675E  -  07 

3. 1676885E  -  09 

1 .5834220E  -  09 

4.1463365E-  11 

2 

1 .81 25565E  -  07 

9.4992060E  -  09 

1.579082  E  -  09 

I.2427645E-  10 

3 

1.8081372E  -  07 

1.5819I51E  -08 

1.5704138E  -09 

2.0674889E  -  10 

4 

1 .6324485E  -  06 

2.8494495E  -  08 

4.2739029E  -  08 

1.1191 609E  -  09 

5 

1 .6304596E  -  06 

8.5448767E  -  08 

4.2621885E  -08 

3.354415  IE  -  09 

6 

1.6264843E  -  06 

1.4229893E  -  07 

4.2387916E-08 

5.58047  E  -  09 

7 

4.52990 18E  -  06 

7.9069730E  -  08 

1  ,f;7742 13E  -  07 

5. 1 780599E  -  09 

8 

4.5243828E  -  06 

2.371 1286E  -07 

1 .97200 13E  -  07 

1.5519987E  -08 

9 

4.51 335 16E  -  06 

3.94867 10E- 07 

1.961 1762E -07 

2.5819375E  -  08 

F 

'4,3 

hS.3 

V 

km"  1 

km"1 

1  1.11  20747E  -  08 

1 . 1090266E  -  08 
1. 1029387E  -  08 
3.0007426E  -  07 
2.9925 178E- 07 
2.9760907E  -  07 
1.3875122E  -06 
1.383709 IE  -  06 
1.3761 134E-06 


2.9 12070 IE  -  10 
8.7282284E-  10 
1.4520463E  -  09 
7.857721  IE -09 
2.3551626E-08 
3.9180977E  -  08 
3.6333286E  -  08 
1.0890027E  -  07 
1.81  16877E  -07 


26 


r-i  c<  r-i 


Time- Independent  Coefficients  of  the  Tide  Potential 


0  0  8.4018456E  -  1 2  4.6140106E  -  1 2  0  0 

l  0  8.368 164 IE  -  1 2  4.5954868E  -  1 2  0  0 

1  1  2.9297877E  -  13  1 .6202500E  -  13  4.3123569E  -  15  2.4544968E  -  1 5 

0  8.3290386E-  12  4.5739464E  -  1 2  0  '  0 

1  2.9177584E  -  13  1.6135948E  -  13  4.2946458E  -  15  2.4444IJ8E-15 

2  2.78405 1 5E  -  1 5  1 .5423222E  -  1 5  8.2132887E  -  17  4.683368  E- 17 

3  0  8.2845357E  -  1 2  4.5494 263E  -  1 2  0  0 

3  I  2.904663 2E  -  13  1 .6063487E  -  13  4.2753632E  -  15  2.4334303E  -  15 

3  2  2.7724443E  -  1 5  1.5358913E  -  15  8. 1790436E  -  17  4.6638388E  -  17 

3  3  1.8513081  E  -  17  1.02S9185E-  17  8.2068535E  -  19  4.6816230E- 

4  3  1.8435 105E  -  17  1 .02 1 5973E  -  1 7  8. 1722863E  -  19  4.6619036E- 
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COMPUTATION  OF  THE  PERTURBING  ACCELERATION 
N*  =  731  +  202  =  933  exactly 
AT  =  .000561214800000  exactly 
dQ  =  28325.500561214800000  exactly 
To  =  .775509940074327173169062286105 
h0  =  .28198651018139069055E  +  05 
s0  =  .37349846089214435790E  +  06 
X  (A/j)  =  2  (h0  -  So) 

=  -  .69059961974801057768E  +  06  deg 
(at*  +  x)  =  ~  .69019706246594063695E  +  06  deg 
=  -  .1204621 1227623637380E  +  05  rad 
cos  (at*  +  x)  =  +  .223888627216 
sin  (at*  +  X)  =  ~  .974614735474 
y,  =  +  316.64861  km 
y>  =  -  6290.36338  km 
y3  =  +  3647. 25332  km 
r  =  7278. 144995  fcm 
p  =  .8763421179 

sin  =  .5011240257 

Uoo  =  54.76683966 
t/,o  =  24.05119117 
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Time-Dependent  Coefficients  of  the  Tide  Potential 


n 

m 

Fnm 

H„m 

0 

0 

-  2.615805043  E-12 

0 

1 

0 

-  2.605292379  E-12 

0 

1 

1 

-  9.231733790  E-14 

-  1.426701083  E- 15 

2 

0 

-  2.593058543  E-12 

0 

2 

1 

-  9.193803464  E-14 

-  1.420837407  E- 15 

2 

2 

-  8.798524745  E-16 

-  2.725617532  E-17 

3 

0 

-  2.579124585  E-12 

0 

3 

1 

-  9.152500570  E-14 

-  1.414451830  E-15 

3 

2 

-  8.761835447  E-16 

-  2.714251175  E-17 

3 

3 

-  5.853884584  E-18 

-  2.725357598  E-19 

Partial  Derivatives  of  Eigenfunctions 

n 

m 

9  u 

dyr 

9  V 
dy, 

0 

0 

— 

.0003273812925 

0 

l 

0 

- 

.0004313144646 

0 

1 

1 

4* 

.006556883732 

+  .0007438816215 

2 

0 

- 

.00009640471307 

0 

2 

1 

+ 

.008605 596947 

+  .001633400725 

2 

2 

4 

.004318488243 

-  .02968388935 

3 

0 

+ 

.0003428662534 

0 

3 

1 

+ 

.001887360387 

+  .001082405391 

3 

2 

+ 

.01195047744 

-  .06493019830 

3 

3 

- 

.1675467571 

-  .03196077254 

n 

m 

a  ,, 

dyi 

a* 

0 

0 

+  .006503572822 

0 

1 

0 

+  .008568250824 

0 

1 

1 

+  .0007438816215 

-  .008183204595 

2 

0 

+  .001915121866 

0 

2 

1 

+  .001633400725 

-  .023760040082 

2 

2 

-  .02585364562 

-  .004125678817 

3 

0 

-  .006811188356 

0 

3 

1 

+  .001082405391 

-  .01956061011 

3 

2 

-  .1057973284 

-  .01400767497 

3 

3 

-  .02763115097 

+  .1246508161 

29 


n 


m 


a 

dy> 


Uy 


dy3 


V.m 


0 

0 

-  .003770875565 

0 

1 

0 

+  .001626320863 

0 

1 

1 

-  .0004313144646 

+  .008568250824 

2 

0 

+  .007577401923 

0 

2 

1 

-  .0001928094261 

+  .003830243731 

2 

2 

+  .03236599779 

+  .003266801449 

3 

0 

+  .005891083242 

0 

3 

1 

+  .00102859876 

-  .02043356501 

3 

2 

+  .04289594100 

-  .004329621566 

3 

3 

+  .02595815241 

-  .1707275267 

Tn  =  -  2.391047870  E-16  km  sec'2 
Ty2  =  -  2.686274153  E-l A  km  sec'2 
Tn  =  -  2.930738682  E-14  km  sec'2 


TXf  =  +  1.467531330  E-14  km  sec~2 
Tn  =  +  2.245096888  E-14  km  sec'1 
TXJ  =  -  2.934580293  E-14  km  sec'1 

Check: 

Tx  =  [Tx]  +  Tx\  +  Tx\]'n  =  3.975659661  E-14  km  sec'2 
Ty  =  [T,\  +  Ty\  +  Ty\\'n  =  3.975659664  E-14  km  sec'2 
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